On the classical-quantum correspondence for the scattering dwell time 
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Using results from the theory of dynamical systems, we derive a general expression for the clas- 
sical average scattering dwell time (r). Remarkably, (r) depends only on a ratio of phase space 
volumes. We further show that, for a wide class of systems, the average classical dwell time is not 
in correspondence with the energy average of the quantum Wigner time delay. 
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I. INTRODUCTION 

The study of the time a quantum collision process takes 
to occur is one of the most interesting chapters in scat- 
tering theory. This problem turns out to be subtle and 
fascinating due to the lack of a Hermitian operator to 
measure the time as a quantum observable. Hence, one 
inevitably has to rely on auxiliary constructions to quan- 
tify the time spent by a scattering process. To that end 
several ingenious strategies have been proposed over the 
past 50 years |1|. In a pioneering work, Eisenbud and 
Wigner Q proposed measuring the scattering delay time 
by recording the peak position of wave packets scattered 
in one-dimension. This simple construction, that just 
invokes the concept of group velocity, already captures 
the deep connection between the energy variations of the 
scattering phase shift and the delay time. In 1960, Smith 
Pi put forward an alternative scheme, applicable to sta- 
tionary scattering processes, where the scattering dwell 
time is associated to the ratio between the probability of 
finding the particle inside the scattering region and the 
flux through its surface. This approach has the advan- 
tage of eliminating the necessity of wave-packets and can 
be easily generalized to multi-channel scattering. As a 
result, the dwell-time tw(,E) is expressed as 
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where the scattering matrix S, that encodes all acces- 
sible information about the scattering process, is taken 
at the energy E. The sums in (^ run over all N open 
asymptotic scattering channels. The time tw{E) is usu- 
ally called Wigner time delay. 

Is t\y{E) in correspondence with the classical dwell 
time for general scattering systems? To answer this ques- 
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tion we approach the problem from the classical side. We 
use the theory of dynamical systems to obtain a remark- 
ably simple and general expression for the classical dwell 
time, revealing its geometric nature. Comparing this re- 
sult with the semiclassical limit for the energy averaged 
quantum dwell-time, we find that the quantum-classical 
correspondence does not hold in general. 

The present analysis does not contradict a previous 
study of ours |j]. There we followed a different path, 
applicable only to chaotic systems, and concluded that 
the classical-quantum correspondence for the dwell time 
holds. Here, approaching the problem in a way that is 
insensitive to the details of the dynamics, we vastly ex- 
pand U and show that the correspondence fails in the 
more general case of systems with mixed phase space. 

The paper is organized as follows. In Sec. |n] we de- 
rive the central result of this paper, namely, a general 
expression for the classical average time delay in terms 
of the system phase space volume. It is key to our anal- 
ysis the formulation of the scattering process as the first 
return of a measure-preserving map, allowing us to ben- 
efit from well-known results of ergodic theory. In Sec. lIIII 
we discuss the semiclassical limit of the energy-averaged 
Wigner time delay. We conclude presenting, in Sec. IIVI 
a comparison between the classical and quantum dwell 
time. We show that, in general, these two quantities do 
not coincide. 



II. THE CLASSICAL DWELL TIME 

Poincare sections are extremely useful tools for the 
analysis of phase space structures in bounded low- 
dimensional Hamiltonian systems: These surfaces of sec- 
tion allow us to reduce the continuous time evolution of 
dynamical systems to discrete mappings, much simpler 
to work with. 

Surfaces of section are essential for the proper defini- 
tion of a scattering problem. Consider the scattering of 
a particle by a potential. The description of the scat- 
tering process requires two control surfaces for detecting 
the state of the particle before and after the scattering 
event. The description of all possible scattering processes 



demands the control surfaces to be chosen so as to en- 
close the scatterer completely. In this case we can con- 
sider just one surface for registering both the states of 
incoming and outgoing particles. 

Let us illustrate these concepts by discussing a generic 
scattering process in three-dimensions. We choose a 
spherical control surface enclosing the region where the 
potential is non-negligible. A point on the associated 
Poincare surface E has coordinates (q, P||), where q rep- 
resents a position on the sphere and p|| the conjugate 
(angular) momentum. An incoming state is completely 
specified by giving its coordinates on S together with 
the condition that the momentum normal to the sphere, 
p j_ , must point inwards (the modulus of p^ is fixed by 
energy conservation). The incoming state then evolves 
inside the scattering region, along a trajectory given by 
Hamilton equations. It eventually intersects S again at 
the exit point (q', pj,) and escapes. Hence, any scattering 
process can be essentially viewed as the first return map 
of E [3, 
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As a consequence of the Poincare-Cartan theorem, this 
map is volume-preserving [f|. 

The structure of the classical scattering problem has 
a clear quantum mechanical counterpart. The quantum 
analogue of the classical Poincare surface is the Hilbert 
space TLy. associated to S. The quantum scattering S- 
matrix is a linear operator of Tis, mapping incoming 
states into outgoings ones. The Poincare map Q is the 
classical limit of S. Conversely, the scattering matrix S 
can be thought of as the quantization of 5 |7| • The uni- 
tarity of S is the quantum counterpart of the classical 
volume conservation [aQ- 

This parallel between classical and quantum scatter- 
ing processes serves to facilitate the determination of 
some quantum-classical correspondences. For instance, 
and very useful for what follows, it becomes clear that 
the classical analogue of an average over "channels" (a 
complete basis set of Tis ) is an average over E weighted 
by its Liouville measure. 

Let us now discuss in detail a very simple scattering 
system: a two-dimensional billiard with an attached pipe. 
The case of a smooth cavity with several (smooth) pipes 
in two or three dimensions, or even, the scattering of 
asymptotically free particles by a smooth potential, are 
conceptually equivalent to the two-dimensional billiard 
with a single pipe, and will be discussed later. 

The physical process we analyze is: A classical parti- 
cle propagates along the pipe and eventually arrives at 
the billiard, where it elastically bounces n times at the 
walls before escaping (see Figure P). A Poincare section 
S, transverse to the pipe, separates the scattering re- 
gion (interior, billiard region, interaction region) from the 
asymptotic region (exterior, pipe). We ask the average 
number of times (n) a particle bounces before escaping, 
or, what is the average dwell time (r) of a particle inside 
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FIG. 1: (top) Billiard with attached pipe. The Poincare sec- 
tion E (dashed line) defines an auxiliary closed billiard, (bot- 
tom) Boundary phase space F of the scattering billiard. The 
phase space coordinates are q, the position along the bound- 
ary, and p, its conjugate momentum. The shaded rectangle 
C corresponds to the closure E. 



the billiard. As already mentioned, the appropriate mea- 
sure for averaging gives equal weights to all points on E 
having the same energy E. 

In what follows we show that the answers to these ques- 
tions are given by very simple ratios between phase space 
volumes. Then we argue that our results are also appli- 
cable to more general geometries. 



A. BirkhofF maps 

Let us consider the Birkhoff section taken along the 
billiard walls, see Fig. ^ The coordinates {q,p), where 
q is the particle position on the billiard boundary and 
p is its conjugate momentum, entirely characterize the 
particle phase space F. The dynamics inside the scatter- 
ing region is given by the Birkhoff (or boundary) map 
T that propagates a particle between successive bounces, 
i.e., from a phase space point {q,p) G F to the one where 
the next bounce takes place. 

Now we "close" the billiard by adding a straight seg- 



merit normal to the pipe axis Q- The Poincare section 
associated to this segment, S, closes the Birkhoff section 
r. Thus, the scattering process can be identified with 
the first-recurrence map to S, now considered as a part 
of F; the dwell time becomes the first-return time to S. 
Figure ^ shows the boundary phase space F, namely, a 
rectangle of length equal to the perimeter of the closed 
billiard and height 2pmax, with p^^x = 2raE, where m 
is the particle mass and E its energy. The shaded ver- 
tical strip corresponds to the closure S, and is denoted 
by C. The inclusion of more pipes to the billiard is ac- 
counted for by adding the corresponding (disjoint) verti- 
cal stripes. This construction can be easily extended to 
higher dimensions. 

We recall that n is the return time measured in units 
of bounces against the billiard walls. Its average is 
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where C„ C C is the subset of initial conditions that first 
return to C after n iterations of the boundary map T 
and /i refers to the volume measure in F. Using measure- 
preservation arguments, it is not difficult to show that 
Er=i'^M(a) = Mdjr^i^"^) [iQ|- Hence, Eq. ® be- 
comes 
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which now expresses the time delay as a quotient of two 
measures: The denominator is the measure of the clo- 
sure; the numerator represents the measure of the inner 
phase space that is explored by the ensemble of scattering 
trajectories. For an ergodic dynamics the set IJtiLi ^"C 
clearly coincides with the full phase space F. Remark- 
ably, even nonergodic Birkhoff's maps very often satisfy 
the weak ergodicity condition 
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For instance, it is simple to verify that the circle bil- 
liard IMIlll) an archetype of integrable dynamics, satis- 
fies Eq. © for any straight closure S. 

As a result, for weakly ergodic billiards we find that 
(n) = /x(F)//Lt(C). The weak-ergodicity condition is not 
verified by systems containing stable islands that cannot 
be reached from the outside, as for instance, the cosine- 
shaped billiard n2\ . In these cases we replace F by F', the 
phase space that is effectively explored by the scattering 
orbits to write 
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The expression above shows that for both ergodic and 
non-ergodic systems (n) is finite. Thus, the probability 
of first returning after n iterations. 
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must decay faster than 1/n?. In some cases numerical 
simulations may suggest a divergent average return time. 
Note, however, that the true asymptotic decay may settle 
only after very long times |l3| | . 



B. Continuous time 

The real, continuous, time-delay problem is addressed 
in analogy to the simple one presented above. To make 
a link between continuous dynamics and maps we invoke 
the stroboscopic map Ta*, i.e., a discretization of the 
continuous evolution into time steps of length At. The 
continuum limit is obtained by making At — > 0. The map 
T^t acts on the full phase space of the scattering system, 
namely, a four-dimensional space for a planar billiard. 

In order to adapt Eq. © to the present context, we 
note that the set of incoming states C has zero measure 
when thought of as a subset of the full phase space: It has 
to be substituted by a properly defined set having finite 
measure, which we call C. The simplest way of choosing 
C is by letting C acquire two extra dimensions: in the 
direction normal to the energy surface, and in the direc- 
tion parallel to the phase space fiow. The corresponding 
additional canonical coordinates are the energy E and 
the time t measured along the trajectories starting at 
the section E. The variables E and t, together with the 
coordinates on E, form a local canonical set. When C 
grows in "thicknesses" by AE and Ai, we have the sim- 
ple relation between measures 



fi{C) = ^i{C) AE At 
The dwell time is then given by [lj| 



(r) = lim At ^ffl , 
where 
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The quantity F' represents the inner phase space for the 
continuous dynamics that can be accessed from outside. 
It has an energy thickness AE. 

By construction, the set C has the important property 
that all its points enter the scattering region after one 
time step Af, namely. 
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This avoids the problem of having to subtract spuri- 
ous contributions to the dwell time arising from non- 
scattering orbits [ij. 

Let us now express Eq. (Q in terms of more appeal- 
ing quantities. Define 17 to be the phase space volume 
contained by the energy shell E within the scattering 



region (as before, primes will indicate "accessible from 
outside"). Then 
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We recall that ij,{C) is the phase space volume contained 
by the energy shell E within the section S. We now 
switch to a more standard notation and, from now on, 
we call it n^. Gathering everything and substituting 
into Eq. (O we arrive at 
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This remarkable formula is exact and holds irrespective 
of the dynamics being chaotic, regular or mixed. After 
making the proper identifications, Eq. p3|l can also be 
applied to billiard systems in three (or higher) dimen- 
sions. 

For the sake of illustration, let us, for instance, use 
Eq. (|13|l to calculate the mean time between collisions 
for a closed billiard. In this case, we have to choose S 
as the phase space corresponding to the full boundary C 
The weak ergodicity condition is obviously satisfied and 
the average bounce time reads t-q — {dVt/dE)/Vlc for any 
billiard, truly ergodic or not. For the two-dimensional 
case, Q,c = 2mvC and il = np^A — 2nmEA {v = ve- 
locity, C = perimeter, A = area). Then we arrive at the 
well-known result Tb = irA/vC, as heuristically shown in 
Ref. 1^ and proven in Ref . jlQ] . 



C. Smooth systems 

The extension of our findings to smooth systems is im- 
mediate. The boundary that so far defined the system 
billiard-plus-pipe now is thought as the level curve of a 
smooth potential. The motion in the waveguide is free 
in the longitudinal direction (x). In the directions trans- 
verse to the waveguide (y) the dynamics is governed by 
a smooth Hamiltonian H±(y,py). 

The analysis of Sec. IIIBI applies equally well to this 
case Thus, the formula for the dwell time is also Eq. (|13|l . 
with the following definitions, f^s is the measure of the 
phase space in the Poincare section lying inside the en- 
ergy shell H± ~ E, 
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fl is the volume of the inner phase space with energy less 
than E. Assuming that the scatterer lies in the region 
X > 0, we have 
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The case of a particle scattered off a smooth potential 
in three dimensions can be accounted for by enclosing the 



scatterer with a large enough spherical shell (the Poincare 
section E). Then one defines the delay time as the (aver- 
age) return time to E minus the return time when there 
is no potential. Both return times are special instances 
of Eq. H13I) . the free-fiight time being just the average 
bounce time of a spherical billiard. 



III. AVERAGE WIGNER TIME-DELAY 

The Wigner time delay tw{E), given by Eq. QJ, fiuc- 
tuates as a function of the energy. Large time delays are 
due to resonant scattering, whereas off-resonance scatter- 
ing corresponds to direct processes that spend short times 
in the interaction region. This picture becomes particu- 
larly clear in the regime of isolated resonances: Long time 
delays occur at narrow energy windows around each res- 
onance, in the remaining energy interval scattering pro- 
cesses are fast (direct). The important energy scale that 
emerges from this picture is the mean resonance spacing. 
When the resonances are overlapping, the separation of 
time scales is less clear, and fiuctuations are much smaller 

By averaging the Wigner time delay over an energy 
window A_E containing many resonances, fast and slow 
processes concur to give a very simple expression: 
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where p{E) is the mean resonance density (due to the 
scattering region). Hence, {tw) basically just counts the 
number of resonances within the energy interval Ai5. 
Equation (|16|l can be derived in various ways, for in- 
stance, using the S-matrix pole structure |17; |. 

In order to relate {tw) with the classical results, we 
take the semiclassical limit of Eq. (|16|) . We first use 
the Weyl formula to express the mean resonance density 
p{E). For that purpose, we consider the corresponding 
closed system (scattering region closed by E), to write 
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where d is the dimension of the system. The wide tilde 
is used to indicate that the semiclassical limit was taken. 
By the same token, the number of states in the pipes is 
given by 



We then arrive at 
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Remarkably, as in the classical case, the average Wigner 
time delay is a purely geometric quantity, and does not 
capture dynamical features. 



IV. CONCLUSIONS 

The most striking result of our semiclassical analysis 
in that the Wigner time delay of Eq. (|19|l is not in corre- 
spondence with the classical dwell time of Eq. (|13|l . The 
correspondence holds only in the case of weak ergodicity, 
where the phase space volume Vl' equals f2. Both quanti- 
ties are different in the more general situation of a mixed 
phase space. 

This result can be interpreted as follows: In general, 
mixed systems have phase space domains in the interac- 
tion region which are not classically accessible from the 
outside. These regions, if larger than h'^ will support 
quantum states. Such states correspond to resonances, 
that can be very thin, depending on the height of the 
dynamical tunneling barriers. As we showed, they con- 
tribute to (tw) with the same weight as other quantum 
states predominantly localized in classically accessible re- 



gions. 

In a broader picture, we speculate that the lack of 
classical-quantum correspondence for the dwell time is 
another manifestation of the non commutativity between 
the long time limit {t —>■ oo) and the semiclassical limit 
{h -^ 0). Tunneling into (or out from) localized states 
at islands of the mixed phase space takes a very large 
time scale to occur, and is absent in the classical limit of 
h = 0. 

We conclude by stressing that our results are rigorous, 
and do not depend on the here presented interpretations. 
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